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Abstract. We consider a random graph model evolving in discrete time- 
steps that is based on 3-interactions among vertices. Triangles, edges and 
vertices have different weights; objects with larger weight are more likely to 
participate in future interactions. We prove the scale free property of the 
model by exploring the asymptotic behaviour of the weight distribution. We 
also find the asympotics of the weight of a fixed vertex. 

1. Introduction 

We consider a random graph model evolving in discrete time-steps. The most 
important feature of the model is the presence of 3-interactions among vertices. 
In our graph process vertices, edges and triangles get nonnegative, integer valued 
random weights which grow with time. The dynamics is driven by these weights. 

The model we are going to deal with resembles those in p], [3], and [4], but 
there are essential differences. In [3] there is no interaction between more than 
two vertices, and the weight of a vertex is simply equal to its degree. In [4] and 
PQ interactions among groups of vertices do appear, but with completely different 
dynamics. 

Our goal is to prove that the ratio of vertices of weight w tends to some positive 
constant x w almost surely, as the number of steps goes to infinity. We will give a 
recursion for x w , from where it will be easy to see the polynomial decay of x w as 
w — > oo. This is the so-called scale free property [2]. We will also determine the 
asymptotics of the weight of any fixed vertex. In the proofs martingale methods 
from [4] are used. 

2. The model 

We start with a single triangle. This has initial weight 1, and all its three edges 
have weight 1. Later on, we will add vertices and edges to the graph randomly. 
Vertices, edges, and triangles will have nonnegative integer-valued weights, which 
increase according to the random evolution of the graph. 



Date: 1 December 2011. 

2010 Mathematics Subject Classification. 05C80, 60G42. 

Key words and phrases. Martingales, random graphs, preferential attachment, scale free 
property. 

The European Union and the European Social Fund have provided financial support to the 
project under the grant agreement no. TAMOP 4.2.1./B-09/KMR-2010-0003. 



2 



Agnes backhausz and tamas f. mori 



The graph is evolving in discrete time-steps. The sum of the weights of triangles 
will be increased by 1 at each step, while the total weight of edges will be increased 
by 3 step by step. 

At each step either a new vertex is added, which then interacts with two already 
existing vertices, or 3 old vertices interact. This has to be decided at the beginning 
of the step, independently of the past. The probability that a new vertex is born 
is p at every step; this is a parameter of the model. We will need < p < 1. 

Assume that in the nth step a new vertex is added to the graph. We choose two 
of the old vertices randomly; they will interact with the new vertex. With proba- 
bility r, independently of the past, the choice is done according to the "preferential 
attachment" rule, and with probability 1 — r it is done "uniformly" . r is a fixed 
parameter of the model. More precisely, in the case of "preferential attachment" 
we choose the endpoints of an already existing edge having weight w with proba- 
bility gj^. Note that the sum of the edge weights is equal to 3n at this moment. In 
the case of "uniform" selection two vertices are chosen with each pair having equal 
probability to be selected; that is, we perform sampling without replacement. This 
allows us to generate edges between old vertices that were not connected before. 

Then the new vertex interacts with the two selected vertices. This means that 
the triangle they form comes to existence with initial weight 1. The two new edges 
connecting the new vertex to the other two get weight 1 each. We connect the old 
vertices if they are not connected yet with an edge of weight 1. This may only 
happen with uniform selection. If the two old vertices are already connected, then 
we increase the weight of that edge by 1. To put it in another way, we increase 
the weights of all three edges of the 3- interaction by 1. This is the end of the step 
where a new vertex is generated. 

With probability 1 — r, 3 of the old vertices will interact. With probability q they 
will be chosen according to the "preferential attachment" rule, and with probability 
(1 — q) they will be chosen "uniformly" . q is the third parameter of the model. This 
choice is also independent of the past. 

In the "preferential attachment" case we choose an already existing triangle of 
weight w with probability proportional to its weight, that is, with probability ^. 
Note that there may exist triangles with zero weight in the graph. A triangle has 
positive weight if and only if it has already appeared in a 3-inter action before. 

On the other hand, in the "uniform" case three distinct vertices are chosen such 
that every triplet has the same probability to be selected. This is again sampling 
without replacement from all the existing vertices. 

In both cases, having selected the three vertices to interact, we draw the edges 
of the triangle that are not present yet. Then the weight of the triangle is increased 
by 1, as well as the weights of the three sides of the triangle. That is, the initial 
weight of a newly generated edge is 1, while the old ones' weights are increased by 
1. 

Now we define the weights of vertices. The weight of a vertex is the sum of the 
weights of the triangles that contain it. Note that this is just the half of the sum of 
weights of edges from it, because whenever a vertex takes part in a 3-interaction, 
the first sum is increased by 1, and the latter one is increased by 2. 

Denote by T n the cr-field generated by the first n steps, and by V n the number 
of vertices after the nth step. Thus Vq = 3. Since we decide independently at each 
step whether a new vertex is born, by the strong law of large numbers we obtain 
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that 

(1) V n = pn + o(n 1/2+e ) a.s. 

for all e > 0. 

Throughout this paper, for two sequences (a n ) , (b n ) of nonnegative numbers, 
a n ~ b n means that b n > except finitely many terms, and a n /b n — > 1 as n — > oo. 

3. Asymptotic weight distribution 

We are interested in the distribution of weights of vertices. As we mentioned 
before, this is the half of the degree of a vertex counted with multiplicity. 

Scale-free property often emerges in models where preferential attachment rules 
are applied. Therefore, throughout the paper we suppose that the parameters do 
not exclude preferential attachment; that is, either r > 0, or q > and p < 1. 

X [n, w] denotes the number of vertices of weight w after n steps. Our goal is to 
examine the asymptotic behaviour of ; more precisely, to prove that the ratio 

of vertices of weight w is convergent almost surely. The limits are deterministic 
constants, which form a polynomially decaying sequence asra-y oo. We may refer 
to this fact as the scale free property of the model, following the terminology of 
Albert and Barabasi [2]. We also compute the characteristic exponent. 

Theorem 3.1. For w — 1, 2, . . . we have 

X[n, w] 

almost surely, as n — > oo. The limits positive constants satisfying the 

following recursion. 

_ 1 _ a(w - 1) + 

a + p + 1 aw + p + 1 



where 



Moreover, 



l P r+(l-p)q>0, P=-[2p(l-r)+3(l-p)(l-q)]. 

6 p 



x w ~ Cw V a I , 
as w — > oo, with some positive constant C . 

Finally we remark that (xd) is a probability distribution, its sum is equal to 1. 

Proof. First we compute the probability that a given vertex of actual weight w 
takes part in the 3-interaction at the nth step. 

If at the nth step a new vertex is generated and we follow the "preferential 
attachment" rule, then this probability is equal to since the total weight of 
edges is 3n, and the sum of weights of edges from the given vertex is just the 
double of its weight w. 

On the other hand, there are (^"j -1 ) pairs of vertices, and every vertex is con- 
tained in V n -i — 1 of them. Hence the probability of being chosen at uniform 
selection is tt^—. 

V n —i 

Now let us examine the case when old vertices interact. A vertex of weight w 
is contained in triangles of total weight w by definition, while the total sum of 
triangles is equal to n after n — 1 steps. Hence the probability of being chosen 
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is — by the "preferential attachment" rule. With "uniform selection" it is clearly 

( y "V i )/(V) = vfc- 



Putting this together we get that the probability that a vertex of weight w takes 
part in the interaction of step n is given by the following. 



(2) P 



2w M ,2 
r-=- + (l-r) 
Sn 



V n -i 



+ (l-p) 



w , ,3 

9- + (l-?hr- 

n V n - 



aw 



dp 



n V n - 



Now we determine the conditional expectation of X[n, w] with respect to T n -\- 
The weights can change at most by 1 . After n— 1 steps we have X [n — l,w] vertices 
of weight w. Each of them increases its weight with the probability given above; 
while vertices of weight w — 1 will count if they take part in the 3-interaction at 
step n. Using the additive property of expectation we obtain that for n > 1, w > 1 
the following holds. 



(3) E(X[n,w] J"„_i) = X[n-l,w] -X[n- l,w 



X[n - l,w- 1] 



n V n -i 
a(w — 1) /3p 



= X [n — 1 , w] 



1 - 



aw j3p 
n V n -i 



X[n - l,w- 1] 



a(w — 1) j3p 



+ p5 w ,i 



- p5. 



where X[0,w — 1] is meant to be zero. The last term is only present for w = 1, 
because the weight of the new vertex is 1. 

We define the following normalizing constants. 



n-l 



n ' 



aw (3p 



V; 



i-1 



n > 1, w > 1. 



By equation ((T|), with any positive e less than ^ we have 

\ 



logc[n,w] = ^ - log 1 



aw 
i 



i + o[i 



l/2+ S \ 



n-l 



E ( ? + t + °( r3/2+£ ) ) = + /3) E - + o(D 



almost surely, where the error term converges as n —> oo. This implies that 

(4) c[n,w] ~ a u ,n Qtu+ ' 3 a.s. 

as n — ^ oo, where positive random variable. 

Introduce Z[n,w] = c[n,w]X[n,w], n > 1, w > 1. From equation ([3]) it is 
clear that (Z[n, w], J- n ) is a nonnegative submartingale for every positive integer 
u). Consider the Doob-Meyer decomposition Z[n,w] = M[n,w\ + A[n,ui], where 
M[n,w] is a martingale and A[n, w] is a predictable increasing process. Based on 
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equation ([3|) we have 

n 

(5) A[n, w] = EZ[1, to] + 5^(E(Z[i, to] | - Z[t - 1, to]) 



i=2 



EZ[1,H + S [X[i-l,w-l] ( a(w , + 



i=2 



i-1 



We will also need a bound on the variation of the martingale part. By using 
equation (QJ we obtain that 



(6) B[n,w) =J2v&r(Z[i,w}\J r i - 1 ) = ^ c[i, w] 2 Vav(X[i, to] \ Fi-i) 

i=2 i=2 
n 

= ^c[i,w\ 2 Vai(X[i,w] -X[i- l,to] | 

i=2 
n 

< ^c[i,<io] 2 E((X[i,tu] - X[i - l^to]) 5 



i=2 



<9^c[z, W ] 2 = 0(n 2 ^+^ +1 ). 



First we used the facts that c[i,w] is measurable with respect to T n -\i and, since 
there is exactly one 3-interaction at each step, the change X of is less than or 
equal to 3. Note that B[n,w] is just the increasing process in the Doob-Meyer 
decomposition of M[n, w] 2 . 

The proof continues by induction on w. For to = lwe obtain that 



(7) A[n,l] ~p^c[i,l] -p^aii Q+/3 ~p : — - ■ n 



a + f3 + l 



almost surely, as n — » oo. 

On the other hand, B[n, 1] = O (n 2 ( Q+/3 ) +1 ) , hence, by applying Proposition 
VII-2-4 of Neveu [5] to M[n, w] we get that 

M[n,w] =o(B[n, 1] 1/2 log 1]) = o(A[n, 1]) 
(see Section 6 of [4] for more details of this argument). Finally we obtain that 

Z[n, 1] ~ A[n, 1] a.s. 

as n — > oo. Using the asymptotics of c [n, 1] and A [n, 1], that is, equations Q and 
([7]), then equation ((T|) and the definition of Z [n, 1], we get that 



X[n,l] _ Z[n,l] _ a + p + l F 1 



V n c[n,l]V n a x n a+ P-pn a + (3 + 1 

almost surely, as n —> oo. 

Thus the theorem holds for w = 1 with a?i 



a+/3+l' 

Suppose that the statement of Theorem 13 . 1 1 holds for to — 1 for some fixed w > 2; 
that is, the ratio of vertices of weight to — 1 converges to some constant x w -\. By 
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using this fact we can compute the asymptotics of A[n, w}. From ([5]) we have 



A[n, w] ~ c[i, w]X[i — 1, w — 1 

«=2 
n 



a(w — 1) ftp 



i Vi-i 
a(w — 1) ftp 



~ o^-i j^p (a(u> - 1) + /3) i QUJ+ ' 3 

z=2 

_ Q^-ipJaCtP - 1) + ft) nau , +/ j+i 
aw + /3 + 1 " 

almost surely, as n — > oo. Here we also used that a and ft are both nonnegative, 
which is clear from their definition. 

From inequality ^ we know that B[n,w] — O [n 2 ^ aw+ ^ +1 ) , thus Proposition 
VH-2-4 of [5] can be applied again to conclude that 

M[n, w] = o (B[n, w} 1/2 logB[n,w]j = o(A[n, id]). 

We end up with 

a w x w -ip(a{w -l)+ft) nC ™ +/3+ i 

^ r , aw + (3 + 1 " a(w-l)+ft 
X [n, w\ ~ 11 -5 = — — np 

L J a w n^+l 3 aw + ft + 1 

almost surely, as n — > oo. Hence 

X[n,w] a[w-l)+ft 

hm — = x w - x ,0,1 a - s - 

n->oo v n aw + ft + 1 

Thus the induction step is complete: x w exists, and it is positive and finite. 
Furthermore, we have a recursion for x w , from where 

(8) x - Jx fr «0--i) + ^ - 1 ff.i±£. 

7 = 1 J 



J=2 



with some positive constant C, as id tends to infinity. The proof of the theorem is 
complete. □ 



( 1 + e+k) T (w + £) ( i\ 



ar(i + ^)r( w + ^±i + i 



Remark. From ((SJ it follows that 



n7 =1 Vi+/3) 

2Af-> Vw — ll 



where 

aj + ft 



n 



A a j + /? + 1 

Zience 

00 

^ =yo = i. 



UJ — 1 
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4. The weight of a fixed vertex 



In this section our goal is to determine the asymptotics of the weight of a fixed 
vertex. Since the weight of a vertex is just the half of its degree when edges are 
counted with multiplicity, we could reformulate our result to obtain the asymptotics 
of the degree. 

It is clear that the weights of the vertices of the starting triangle are interchange- 
able, therefore it is not necessary to deal with all the three. Let only one of them 
be labelled by 0, the other two will remain unlabelled. Moreover, let the further 
vertices get labels 1, 2, etc, in the order they are added to the graph. Let W [n,j] 
be the weight of vertex j after step n, provided it exists. Otherwise let W[n, j] be 
equal to zero. Obviously, vertex j cannot exist before step j. Let I[n, j] denote the 
indicator of the event {W[n, j] > 1}. 

Let us introduce the sequences 

h 



d n = ftp 



i=l i=l ' " 

with a, ft defined in Theorem l3.ll Note that b n is deterministic, while d n is random, 
but J-" n _i-measurable. 

Lemma 4.1. Let j and k be fixed integers, < j < k, and let Z[n,j] = 
b n W[n,j] — d n . Then (Z[n,j]I[k,j], J-" n ) is a martingale for n > k. 

Proof. According to equation @, the probability that it gets new edges at step 



1 is equal to 



aW[n,j) ftp 



, provided that vertex j already exists, which surely 



n + 1 V n 

holds if the indicator I[k,j] differs from 0. This implies that 

aW[n,j] 
n + 1 



a 



E(l[k,j]W[n + l,j] | F n )=I[k,j]W[n,j] + I[k,j] 

= I[k,j]W[n,j] ll- 
Multiplying both sides by b n +±, we get by definition that 

E(b n+1 W[n + l,j]I[k,j] | F n ) = I[k,j] ( b n W[n,j] + b n . 



ftp 



I%3 



ftp 



ftp 



Vn, 

= I[k,j](b n W[n,j] -d n + d n+ i), 

which completes the proof of the lemma, since d n +\ is J^-measurable. □ 

Theorem 4.1. Fix j > 0. Then W[n, j] ~ Cj nC * almost surely «sn-> oo, where 
Q is a positive random variable. 

Proof. First we show that this holds with a nonnegative Q. 

Almost surely on the event that vertex j exists after step n we have 

¥(W[n + l,j] = W[n,j] + 1 | F n ) > . 

Using the Levy-type generalization of Borel-Cantelli-lemma [5, VII-2-6] we get that 
VF[n,j] — > oo with probability 1. 

From the definition of b n it easily follows that 

r(n + i)r(i + a) 



(9) 



b n = 



T(n + 1 + a) 



r(l + a)n _a , 



s 
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asn-> oo. Hence, by using equation (JlJ and the positivity of a, we get 

n , n — 1 

rf « = & E if- = p ^ + a ) E r C 1 + ■ 



1=1 



Thus d„ is almost surely convergent as rt — > oo, hence the martingale of Lemma 
4 . 1 1 is bounded from below. This martingale has bounded differences, for 

Z[n + 1, j] - Z[n,j] < b n (W[n + 1, j] - W[n,j}) < b n < 1, 



and 



Z[n,j] - Z[n + l,j] < (b n - b n+1 )W[n,j] + (d 



dp 



< b n +ia + b n+1 ^- < (« 



3 



By Proposition VII-3-9 of [5] such a martingale either converges or oscillates be- 
tween — oo and +oo, but now the latter is excluded, hence it must converge almost 
surely. 

Going further, we get that is convergent almost everywhere on the 

event that the weight of vertex j is greater than 1 after step k. Since that weight 
tends to infinity, the limit as k —¥ oo of this increasing sequence of events has 
probability 1. Thus b n W[n,j] is almost surely convergent, and by equation (j9|) we 
get that the statement of the theorem holds with a nonnegative Q. 

Now we have to prove that Q is positive with probability 1. 

In what follows, if the indicator in the numerator is zero, let us define the frac- 
tions to be zero. Similarly to the previous lemma, for n > k we can write 



E 



I[k,j] 



W[n + l,j]-l 



aW[n,j] | (3p\ I[k,j 

Vt. 



1 



1 



It is clear that 

aW[n,j] (3p 
n+1 V„ 

hence we get that 



I[k,j] 



I[k,j] 

W[n,j] W[n,j]-1 



W[n,j] 

aW\nJ] 
n+1 



< 



v„ 



I[k,j] 



W[n,j}-1 
al[k,j] 



(n+l)(W[n,.7]-l)' 



I[k,j] 



< 



I[k,j] 



W[n+l,j]-l 
From this it follows that 

/ e n I[k,j] 
\W[n,j]- 

is a supermartingale for n > j, where 



W[n,j]-1 



yr( 1 _ ^y 1 = r(l-a)r(n + l) 
J-l \ i) Y(n + l-a) 

4 — 1 V ' 



T(l -a)n a . 



This supermartingale is nonnegative, hence it converges almost surely. Since lim^oo I[k, j] = 
1 holds a.s., we obtain that tttt — — r is also convergent almost surely asn-> oo. 



W[n,j]-1 

This implies that Q > with probability 1, as stated. 



□ 
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